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We consider the coherent control of a quantum bit by the use of short pulses with finite duration Tp . 
By shaping the pulse, we perturbatively decouple the dynamics of the bath from the dynamics of the 
quantum bit during the pulse. Such shaped pulses provide single quantum bit gates robust against 
decoherence which are useful for quantum-information processing. We extend previous results in 
two ways: (i) we treat frequency modulated pulses and (ii) we pass from time-independent baths 
to analytically time-dependent baths. First- and second-order solutions for tt and 7r/2 pulses are 
presented. They are useful in experiments where amplitude modulation is difficult to realize. 



PACS numbers: 03.67.Pp, 82.56.Dj, 76.60.Lz, 03.65.Yz 

I. INTRODUCTION 

The occurrence of decoherence in quantum systems 
is one of the main difficulties to be overcome in mod- 
ern experiments. Especially quantum information pro- 
cessing (QIP) requires the quantum system to remain in 
fixed phase relations between the application of quantum 
gates. Otherwise one can not benefit from the quantum 
parallelism making QIP such a powerful technique. 

A generic example for a quantum bit (qubit) is a spin 
S — 1/2. One of the first implementations of a quantum 
algorithm was realized with nuclear magnetic resonance 
(NMR) [T]. The qubits are encoded in the nuclear spin 
degrees of freedom. We also use the spin language to 
describe the states and dynamics of the qubit. The state 
t is identified with the logical |1) and the state I with 
the logical |0). 

The loss of coherence is induced by the coupling of 
the spin to its environment, the so-called bath. One way 
to suppress this coupling is the application of suitable 
control pulses first introduced by Hahn in 1950 2 for 
NMR experiments. This idea led to the development 
of sequences for control pulses ranging from the Carr- 
Purcell-Meiboom-Gill (CPMG) cycle [1 H] to more and 
more complex control schemes ^ . In QIP, this approach 
is known under the name of dynamic decoupling (DD) 
[BHH] . Theoretically, dynamic decoupling can be achieved 
to infinite order in the duration T of the total pulse se- 
quence OdO]. One particularly efficient way to deal with 
pure dephasing decoherence is the use of theoretically op- 
timized DD (Uhrig DD) pTMl3| . It has been successfully 
implemented experimentally |14L I15j . 

The pulses considered in theoretical studies of DD 
schemes are mostly ideal in the sense that they have an 
infinite amplitude and that they act instantaneously in 
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time corresponding to Dirac 5 pulses. Of course, this 
property can not be realized in experiments. If the fi- 
nite pulse duration is taken into account in calculations, 
it turns out to be a nuisance in most cases (see, for in- 
stance, Ref s . ITOl and [T5)) . Hence, there is an abundant lit- 
erature on pulse shaping and optimization which we can 
mention only partly (for a book see Ref. J^S]). We 

emphasize, however, that a suitably shaped pulse can be 
integrated into a DD sequence such that the high-order 
suppression of decoherence is hardly hampered [331 134) . 

The majority of the existing theoretical studies of pulse 
shaping consider pulses acting on the two-dimensional 
Hilbert space. The goal is to design robust pulses which 
tolerate a maximum of frequency offset or other inaccu- 
racies of the pulse ■ The next stage of com- 
plexity includes random time-dependent classical noise, 
which is still described by classical fields coupled to the 
spin |25j . The maximum stage of complexity considers 
a fully quantum mechanical bath which means that the 
qubit is coupled to a macroscopic quantum mechanical 
system by noncommuting operators [ini [HI E71 - I5T] . It is 
on this level that our present study is situated. We stress 
that a quantum mechanical pulse, which is robust against 
a coupling to its environment, constitutes an appropriate 
single-qubit gate. 

In particular, we extend previous work |30| in two 
ways: (i) We allow for analytically time-dependent bath 
operators, both in the spin-bath coupling and in the 
Hamiltonian of the bath. Such time dependence may, 
for instance, arise from a time-dependent reference frame 
|35) . (ii) We propose frequency-modulated pulses while 
before only amplitude-modulated pulses were studied 
[MIMlllZllMllSnilSI] except in the general no-go theorem 
in Ref. HH] We stress that in the NMR context amplitude 
and phase-modulated pulses have been discussed inten- 
sively [21n23) . But to our knowledge these investigations 
do not comprise quantum mechanical baths nor dynamic 
classical noise. 

Explicitly, we compute continuous solutions for tt and 
7r/2 pulses realized by frequency modulation (see also 
Ref. . The consideration of frequency modulation is 



2 



motivated from experimental situations where the fre- 
quency of a pulse can be controlled more accurately or 
more easily than its amplitude. Thus the present study 
is complementary to preceding ones. 

The paper is organized as follows: In Sec. [llj we give 
an overview of the model under study and motivate our 
ansatz for the time evolution of the whole system. In Sec. 
Ill we derive the perturbative expansion for a generic 
time-dependent bath. We require that the time depen- 
dence is analytical in order to be able to apply a per- 
turbative approach. Then we introduce the frequency- 



modulated ansatz in Sec. IV and specialize the general 
equations for this specific case. The solutions found for 
first- and second-order pulses are discussed in Sec. [V| and 
we finally conclude in Sec. |VI| 



II. MODEL AND ANSATZ 

We consider the general case of a spin coupled to a 
time-dependent bath 



H{t)=H^{t) + <7-A{t), 



(1) 



where H\^{t) denotes the part of the Hamiltonian that 
acts only on the bath. We refer to it as the bath Hamil- 
tonian. The vector of Pauli matrices a acts on the Hilbert 
space of the spin S* = 1/2 while A{t) is a vector of 
bath operators to which the spin is coupled. No spe- 
cial operator structure is assumed for the bath opera- 
tors; i.e., the commutators [Ai{t'), Ai{t)\, [Ai{t'), Aj{t)\ 
and [Ai{t'),H\j{t)] do not need to vanish. 

The model comprises typical cases such as a bosonic 
bath or a spin network. Relevant experimental systems 
comprise the electronic spin in a quantum dot coupled to 
the bath of nuclear isotope spins |36j or the spin of a ni- 
trogen vacancy center in diamond interacting again with 
a bath of nuclear isotope spins [37] • For our purposes, we 
require that the time dependence of the operators H\^(t) 
and A{t) is analytical so that they can be expanded in 
time: 



A{t) = Ao + Alt + A2t^ + ■■■ 



(2a) 
(2b) 



This analyticity is often fulfilled, e.g., in rotating refer- 
ence frames or in the operator interaction picture. For 
fast time-dependencies, however, the above expansion is 
not useful because the derivatives are large. Very fast os- 
cillatory time-dependencies are better treated by average 
Hamiltonian theory. 

Our model includes the common case of a purely de- 
phasing bath, i.e., a spin coupled only along the cr^ direc- 
tion to the bath. This model is justified in experiments 
where the dephasing time T2 is significantly lower than 
the longitudinal relaxation time Ti. This is the case if 
the energetic splitting between the states with cr^ = 1 
and (Tz = —1 is large. 



The coupling strength between the spin and the bath 
is given by A := ||yl(t)|| while the energy of the bath is 
defined to be ojb := ll-H^b(i)||. If these operators are not 
bounded, that means if A and Wb cannot be defined by the 
operator norms, we refer by A and Wb to the generic en- 
ergy scales of the corresponding operators. For instance, 
in a bosonic bath Wb is the upper cutoff of the bosonic 
energy spectrum. The energy scales serve as reference 
values for Tp. That means that we aim at an expansion 
in the dimensionless ratios uj^r-p and Axp. 

Applying the control pulse to the system, the term 



v(t) 



(3) 



is added to the Hamiltonian H{t). Here v{t) is a vector 
of amplitudes describing the controllable time-dependent 
shape of the pulse. The normalized vector v{t)/ \v{t)\ is 
the current axis of rotation at time t while the norm 
\v{t)\ describes the magnitude of the control term which 
determines the velocity of rotation. Without loss of gen- 
erality, we assume that the pulse starts ai t — Q and ends 
at f = Tp. The time evolution between and Tp of the 
combined system reads 



U{t) = T 



exp 



H{t)<it 



(4) 



where T stands for the standard time ordering. 
Our aim is to perturbatively decouple the time evolution 
of the spin from the time evolution of the bath during the 
pulse. This motivates the following ansatz for the time 
evolution of the whole system 



C/(Tp) - t/b(rp)F(rp)C/e(rp), 



(5) 



where 



Ub{t) = Texp {^i Hi,{t')dt'^ , (6a) 

P{t) ■= Texp (^-ia ■ v{t')dt'^ . (6b) 

The unitary operator Uh{Tp) describes the time evolu- 
tion of the bath and P{Tp) the rotation of the spin due 
to the pulse. Note that the ansatz Uh{Tp)P{Tp) does not 
comprise any coupling between spin and bath. It is close 
to the goals of many previous studies aiming at robust 
pulses P"7HS51 m] and it corresponds to the ansatz 
used in previous studies separating the pulse from a class- 
cial [25] or a quantum mechanical dynamics of the bath 
[m [111271 1301 El- We emphasize that an ansatz which 
separates the pulse from the dynamics of the spin plus 
bath system can be shown not to succeed beyond leading 
order [28, 23]. 

Since the spin-bath coupling is not included in 
UhiTp)P{Tp) we introduced the correction unitary opera- 
tor Uc{Tp) in Eq. We want to shape the pulse so that 
the correction term is as close to the identity as possible. 
A perfect decoupling would imply Uc{Tp) = 1. But this is 
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unrealistic to achieve. Hence we pursue the perturbative 
approach to make as many terms of an expansion in Tp 
as possible vanish. Then U{Tp) « Uy,{tp)P{tp) represents 
a valid approximation and one can neglect the spin-bath 
coupling during the pulse. We remark that pulses shaped 
in this way constitute robust single-qubit gates. 



III. DERIVATION 

The derivation of the perturbative conditions for the 
shaped pulses is very similar to the derivation given in 
Ref. 1301 Yet we present a brief outline here in order 
to keep the present article self-contained and because 
we extend the previous derivation to analytically time- 
dependent baths. We start from the pulse Hamiltonian 
in Eq. We describe the time-dependent pulse opera- 
tor as a global rotation about the axis a{t) 



P(t) — exp —ia ■ a{t) 



(7) 



where \a{t)\ = 1. The spin is turned by the angle ip{t) 
at the time t. Every unitary operator acting only on the 
Hilbert space of the spin can be written in the form of 
Eq. ([t]). In particular, a pulse that turns the spin by an 
angle x satisfies 



^(Tp) 



(8) 



We stress the difference between the current axis of rota- 
tion v{t)/ \v{t)\ and the effective axis a{t) describing the 
total rotation of the spin from its position at time to 
its current position at time t. 

By definition the pulse operator fulfills the Schrodinger 
equation 



which implies |29j 



2v{t) = i]/{t)a{t) + a'{t) sva^{t) 
- (1 - cosi/)(t)) [a'{t) X a{t)] 



(9) 



(10) 



This differential equation is solved numerically for the 
frequency-modulated ansatz below. The time evolution 
of the whole system is given by 

idtU{t) = [i/b(0 + S ■ A{t) + U{t). (11) 

Inserting the ansatz ([5| and solving for dtUc{t) yields 

idtUcit) - G(t)C/c(0 (12) 
G{t) = P~\t)U^\t) a ■ Ait) Uy,{t)P{t). (13) 

Thus the unitary correction is determined by a 
Schrodinger equation with G{t) as its time-dependent 



Hamiltonian. The formal solution of Eq. ( 12 ) in terms of 
the standard time ordering operator is 



Uc{t) = T ^exp ^-i^ G(T)dT 



(14) 



Aiming at an expansion of Uc{t) in powers of Tp it is 
convenient to use the Magnus expansion [38^ to express 
the time-ordered exponential 

C/c(rp) = exp (-iTp (g'i) + G(2) + ••.)), (15) 
where G*^'^ — O (xp^^). The first two terms read 

(16a) 



T^G^i) = / dtG{t) 



2 



ti 



C(^)^--l dh dh[G{h),G{h)]. (16b) 
Jo 



Next, we need an expansion of G{t) in powers of time. 
To this end, we consider the representation ([?]) which 
implies 

P^^it)a ■ A(t)P{t) = cos-ipA - sin-0 (a x A^ 

+ (1 - cos ^j) (^a- Aj a -a (17a) 

= ^A(t) (<)•<? (17b) 
= '^n,^j{t)Aj{t)a„ (17c) 



where the time dependencies on the right hand-side of 
Eq. (17a) are omitted to lighten the notation. The vector 



operator nyi(t) (t) is the vector A(t) after a rotation about 
the axis a{t) by the angle —^(t). The corresponding 
rotation matrix Da{—i)) is given by its matrix elements 
nij{t); for their explicit form see Appendix [C[ Due to 
the orthogonality of Da{—ip) the moduli of all its matrix 
elements are bounded by unity. 

Note that there are two different kinds of time depen- 
dence in Eq. (17b). On the one hand, the time depen- 



dence of A{t) becomes weaker and weaker as the pulse 
duration Tp is taken to zero because we assume that A{t) 
is analytical. This is exploited below. On the other hand, 
the time dependence of the ni j{t) scales with Tp, which 
means that hij{s) := nij (sTp) is completely independent 
of Tp because the pulse is completed at t = Tp whatever 
the pulse duration is. 

We proceed by introducing the vector operator A{t) 
and expanding it in powers of t 

A{t) := U-\t)A{t)U^{t) (18a) 
= Ao + it\H^,fl,Ao]+tAi + 0{t^). (18b) 

Here the main differences to the derivation in Ref. [301 
arises. In Ref. [3D1 the term proportional to Ai did 
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not appear because the bath was considered to be time- 
independent. Using the vector operator riy!t(t)(i) from Eq. 
(17b) we rewrite G{t) concisely as 



Git) = p-\t)a ■ Ait)P{t) 



{t)-S. 



(19a) 
(19b) 



This form of G{t) can be expanded in powers of t such 
that the neglected terms are of second order in Tp for 

t € [0, Tp] 

G{t) = Hao ■ S + t{i [fl"b,o, ■ nAo\ + n^i • ct) + O (xp) . 

(20) 

Note that the time dependence stemming from the pulse 
rotation is not expanded because it does not change on 
Tp — > 0. Physically this means that one can expand 
in TpH(t), i.e., in Arp and in WbTp, but not in TpHy{t) 
because the magnitude of Hy{t) has to be increased on 
Tp — > to realize the desired pulse. 



Inserting Eq. ( 20 ) in the terms of the Magnus expan- 



sion (|16p eventually yields 
exp 



■ JQ 
(2) f (2a 

i 

where Ty*-*^ cx t*. Explicitly, one has 



(2c) 



(21a) 
(21b) 

(21c) 



,(2i) 



tnij{t)dt, (22a) 



y^e»jfc»jv(^i)»fc,m(^2 



dti / dt2 
Jo 



= ^ [Aj,Q,Ak.o] r dti f dt2 
inij{ti)ni^kit2) - "ij(i2)"i,fe(ti)) 



(22b) 



(22c) 
and the 



In these equations the anticommutator 
completely antisymmetric Levi-Civita tensor eijk appear. 
The indices k, /, and m take one of the values x, y, or 
z. Most of the Eqs. (21b) and (22) are identical to those 
obtained in Ref. 



Only in rjf"^ does the time depen- 
dence of the bath appear additionally. It is encoded in 
the operators Aj^i which are zero for a time-independent 
bath. 



The second-order equations ( 22 ) do not show a depen- 
dence on pure bath terms of an order higher than -ffb.o- 
This means that the time dependence of the pure bath 
Hamiltonian is irrelevant up to second order. This is 
reasonable because even for a time-independent bath the 



First-order pulses 



FM-l-PI 



FM-1-PI2 



Vo 


3.75146609 


Vo 


4.92892484 


bi 


0.00011442 


bi 


0.00009874 


b2 


-1.09347112 


b2 


-0.94331659 


bs 


0.00012443 


bs 


0.00002530 


b4 


-0.59452572 


bi 


-0.12087663 



Table I. Overview of the pulses satisfying all first-order equa- 
tions (21b I. FM-l-PI denotes the frequency modulated tt 



pulse. FM-1-PI2 denotes the frequency-modulated tv/2 pulse. 
The dimensionless coefficients b„ belong to the ansatz in Eq. 
(33 1. The amplitudes Vo are given in units of 1/rp. With all 
eight digits given the conditions are fulfilled for the tt pulse 
within 10~^° and for the n/2 pulse within 10~^. With only 



two digits they are fulfilled within 10 
tively. 



and 10" 



respec- 



actual bath dynamics induced from i/b appears only in 
second-order conditions. This does not apply to the vec- 
tor operator A{t) which is already relevant in first-order 
pulses. 

In the general case of a completely generic bath, all 
the expressions f?*-"-', a € {1, 2a, 26, 2c}, have to vanish 
in order to fulfill Uc — 1 + O (Tp). The pulse shape de- 
termines the time evolution of the matrix elements Uij. 
Hence, in order to fulfill all the conditions the operator- 



independent integrals in Eqs. (21b) and (22) must disap- 
pear. The resulting 39 scalar equations are identical to 
those obtained in Ref. [301 

This is our first key result. It proves the applicability 
of the previously obtained pulses even in the presence 
of a non-trivial time dependence of the bath which may 
stem from special reference frames or from the interaction 
picture of fast modes. For specific cases, such as the pure 
dephasing model or if [Hi,{t) , A(t)] — the number of 
scalar equations to be fulfilled for Uc — 1 + O (Tp) is 
reduced significantly. Pulses with less complexity can be 
used. This is studied in the sequel. 



IV. FREQUENCY MODULATED ANSATZ 



To solve Ec^. (10) we choose an ansatz for v{t). In this 
paper we focus on a frequency-modulated pulse acting 
only in the axTy-jAane with a fixed amplitude Vq := \v\ 
and the current axis of rotation 



• cos($(<)) \ 
J{t) ={ Vo- sin($(t)) 1 



(23) 



where $(t) is a time-dependent phase. Note the differ- 
ence to pulses with a time-dependent amplitude and a 
fixed axis discussed, for instance, in Refs. fT6 l \T7 \ and[30l 
We focus here on frequency modulation in complement to 
previous work because there may be experimental setups 
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Second-order pulses 



FM-2-PI 



FM-2-PI2 



T / 

Vo 


12.0J432979 


vo 


12.25ol9390 


Ol 


U.ii4(0ioy 


Ol 


i. (30 < io4U 


02 


U.i / 24800 1 


7 

02 


0. (3o299o9 


L 

03 


U.4ozDzozi 


03 


U.Zoz4zozo 


fe4 


-1.14494851 


&4 


-0.24829310 


fes 


-0.20879091 


&5 


-0.07102204 


&6 


0.25378013 


&6 


-0.13192380 


&7 


0.20306835 


&7 


1.07948226 


&8 


-0.16748022 


bs 


0.12220006 


bg 


-0.32052254 




0.04608986 


&10 


0.32586203 


bio 


-0.15365617 



Table II. Overview of the pulses satisfying all first- and 
second-order equations, Eqs. ( |21b|22[ ). FM-2-PI and FM-2- 
PI2 denote the frequency-modulated vr pulse and n/2 pulse, 
respectively. The dimensionless coefficients bn belong to the 
ansatz in Eq. ( 33 1 . The amplitudes Vo are given in units of 



l/rp. With all eight digits given the conditions are fulfilled for 
the TT pulse within lO"'^" and for the 7r/2 pulse within 10~^^. 
With only two digits they are fulfilled within 10"^ and 10"^ 
respectively. 





Minimized second-order 
FM-2-MIN-PI 


pulses 

FM-2-MIN-PI2 


Vo 


10.70711454 


Vij 


8.43541412 


bi 


0.00002087 


fei 


-1.82041507 


b2 


1.38768938 


b2 


-0.35249197 


bs 


-0.00019922 


bs 


0.03054874 


b4 


-0.70668998 


b4 


0.52093576 


fes 


-0.00001588 


fes 


-0.55504440 


be 


0.13773085 


fee 


-0.38815568 


br 


0.00008770 


br 


0.45167361 


bs 


0.68894331 


bs 


-0.19445080 


b9 


-0.00011408 


bo 


-0.16194806 


bio 


-0.69744086 


bio 


-0.28223330 


bi4 


0.46501991 


fel4 


0.04585897 



Table III. Overview of the pulses satisfying all first- and 
second-order equations, Eqs. ( 21b|22 \ with minimized am- 
plitude. FM-2-PI and FM-2-PI2 denote the frequency- 
modulated TV pulse and 7r/2 pulse with minimized amplitude, 
respectively. The dimensionless coefficients fe„ belong to the 
ansatz in Eq. ( 33 1 . The amplitudes Vo are given in units of 



l/rp. With all eight digits given the conditions are fulfilled for 
the TV pulse within 10"^'' and for the tv/2 pulse within 10~^. 
With only two digits they are fulfilled within 10"^ and 10"\ 
respectively. 



where frequency modulation is much easier (or more ac- 
curately) implemented than amplitude modulation. Note 
that the ansatz ( 23 1 assumes that the control pulse can be 
switched on instantaneously. Transients are assumed to 
be sufficiently steep to be taken as jumps. The considera- 
tion of continuous amplitudes and frequency modulation 
is left to future research. 



To point out the relation of the ansatz (23) to the 
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Figure 1. (Color online) First-order vr pulse FM-l-PI. We 
also plot sin<I>(f) oc Vy{t) and cos$(t) oc Vx{t) to illustrate 
the pulse shape in spin space. The left scale refers to $(4) 
and the right scale to sin<l>(t) and cos<I>(t), respectively. The 
Fourier coefficients for this pulse are given in Table |I] 



(24) 



consider a spin with a Larmor frequency ojl in the NMR 
language [55] : 

cr,. 

2 

Of course, this description is not restricted to nuclear 
spins. Any two-level system with an energy splitting can 
be considered. The control field is realized by applying a 
field perpendicular to the cr^-axis rotating with the Lar- 
mor frequency: 

H,[ = Vo {cr^ cos [oji^t - $(i)] - (Ty sin [uji^t - ^{t)]} . 

(25) 

We include a time-dependent phase $(t) to shape the 
pulse. Its derivative dt^{t) is the deviation of the fre- 
quency from the Larmor frequency. In this sense Eq. 



(251 describes a frequency-modulated pulse. Next, H,-{ 
is transformed into the rotating framework in which Hz 
vanishes. Using the unitary time evolution induced by 

Hz 



Urot{t) = exp {i 
we obtain H^ot{t) = u}^^.H^iU^ot, which reads 



i—ta^ 



(26) 



In order to find a{t 
parametrization in Eq, 



Vo 


cos($(i)) \ 




Vo 


sin(<i>(t)) 


(?) 




/ 





t) a 



v{t) ■ 3. (27) 



and ?/;(0 appearing in the 
of the pulse one has to solve 



the difi^erential equation (10). Because a(i) is a unit vec- 
tor, it is convenient to describe it by two angles: (^(i) 
and B{t) 



experimental realization in the laboratory framework we 



sin(6'(i))cos((y9(t)) \ 
a(t) = I sin(6'(i))sin((^(t)) 

cos(6i(t)) ) 



(28) 
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Solving Eq. (10 1 for the time derivatives of ip{t), (p{t), 



and 9{t), we find 
dtip — 2Vo sin 9 [sin $ sin ip + cos $ cos ip] 



dtp = Vo 



dtO = Vo 



cos ^ sin(<i> - (^) 



sin 1^ cos 6' cos($ 



(29a) 



sin |- sin 



(29b) 



cos ^ COS cos($ 



- sin 1^ sin($ — (p) 



(29c) 



The seeming singularities for vanishing angles on the 
right-hand sides of Eqs. (29b) and (29c I have no physical 



reason, but they only result from the choice of spheri- 
cal coordinates and from the chosen parametrization in 
Eq. ([7|. Note that the global axis of rotation a is not 
uniquely defined if is a multiple of 2tt. 

At the very beginning at t = the current axis of 
rotation v and the global one a coincide. The former lies 
by construction in the tTa;CTy-plane. Hence we have the 
initial conditions 



lim^(t) = -, 
lini'0(t) = 0, 
\im(p{t) 

t— i-O 



$(0), 



(30a) 
(30b) 
(30c) 



where the latter two equations represent our deliberate 
choice. Inspecting the limit i -> one additionally finds 



2dtp \t=o^ dtm\t=o' 
dtO |t=o = 0. 



(31a) 
(31b) 



In 



The derivative dt^ follows trivially from Eq. (29a) 
the next section, we provide solutions for this ansatz and 
a specific case of spin-bath coupling. 



V. RESULTS 

We are interested in the experimentally important case 
of a purely dephasing model, i.e., a bath coupled only via 
az to the spin 



A{t) = A{t)ez 



(32) 



Hence the coupling becomes simpler, but the bath dy- 
namics itself is still kept in full generality. Spin flips 
do not occur in this model so that Ti is infinite. But 
decoherence of the T2 type is entirely kept. This as- 
sumption is justified in many experimental realizations. 
Moreover, the simplification of the coupling is advanta- 
geous for pulse shaping because it reduces the number of 
integral conditions derived from Eqs. (21b) and (22) in 



second-order to be fulfilled from 39 to 3 first-order con- 
ditions and 6 second-order conditions which are given 
explicitly in Appendices [X] and |Bj 
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Figure 2. (Color online) First-order n/2 pulse FM-1-PI2. We 
also plot sin<l>(f) oc Vy{t) and cos$(t) oc Vx{t) to illustrate 
the pulse shape in spin space. The left scale refers to $(4) 
and the right scale to sin<l?(t) and cos$(t), respectively. The 
Fourier coefficients for this pulse are given in Table |I] 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

t [units of Tp] 

Figure 3. (Color online) Second-order tt pulse FM-2-PI. We 
also plot sin"l>(f) oc Vy{t) and cos$(t) oc Vx{t) to illustrate 
the pulse shape in spin space. The left scale refers to $(t) 
and the right scale to sin$(t) and cos$(t), respectively. The 
Fourier coefficients for this pulse are given in Table [IT] 



In the following, we present continuous pulses which 
fulfill the first-order conditions (first-order pulses) and 
pulses which fulfill all first- and second-order conditions 



(second-order pulses) for pure dephasing as in Eq. ( 32 1 . 
Thereby, we provide optimized pulses that decouple the 
spin from the bath during the duration of the pulse up 
to 0{t^). 

In order to consider a continuous frequency modulation 
we use the Fourier series ansatz 

<i>(t) = ^62„-i sin {2TTnt/Tp) + b2n [cos {2imt/Tp) - 1] 

n 

(33) 

for ^(t). We consider tt and ti/2 pulses because of their 
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t [units of Tp] t [units of Tp] 



Figure 4. (Color online) Second-order 7r/2 pulse FM-2-PI2. 
We also plot sin$(i) oc Vy{t) and cos<E>(f) oc Vx{t) to illustrate 
the pulse shape in spin space. The left scale refers to 
and the right scale to sin$(t) and cos$(t), respectively. The 
Fourier coefficients for this pulse are given in Table [HI 



Figure 5. (Color online) Minimized second-order tt pulse FM- 
2-MIN-PI; sin4>(t) oc Vy{t) and cos$(t) oc v^{t) show the 
pulse shape in spin space. The left scale refers to $(t) and the 
right scale to sin"l>(t) and cos"I>(f), respectively. The Fourier 
coefficients for this pulse are given in Table [Till 



frequent use in QIP and NMR. Therefore, the pulse has 
to fulfill 



( 35 1 five parameters are necessary to construct first-order 



V'(Tp) = 71" 



tt/2 



(34) 



for TT pulses and 7r/2 pulses, respectively, according to 
Eq. (|8|. The value 0{Tp) is fixed by the fact that the 
final axis of rotation has to be perpendicular to CTz to 
rotate the spin by the full angle ijj{Tp). Thus we require 



^(^p) = o • 



(35) 



pulses. One parameter is the amplitude Vo and the oth- 
ers are the coefficients fo„ in ansatz ( [33| . The character- 
istics of the pulses are reported in Table |l| The pulses are 
plotted in Figs. [T] and [2] Note that the composite and 
continuous amplitude-modulated pulses found in Ref. 1301 
have comparable amplitudes for first-order pulses. 



B. Second-order pulses 



Second-order pulses additionally have to fulfill condi- 



For a given ansatz, Eq. ([33|, the numerical procedure tions (l22]). These equations again simplify for a purely 



to find solutions is straightforward. We solve the differ- 



ential equations (291 using a fourth-order R unge -Kutta 
algorithm. For this solution the conditions (21b), (22|, 
(Isl), and |35|) are evaluated. We search for roots using 



the Powell hybrid method in the GNU scientific library 

Of two different pulses the one with the lower ampli- 
tude Vq is preferable in experiment because less power 
is needed to realize it. For an experimentally realiz- 
able maximum amplitude this implies that the theoretical 
pulse with lower amplitude can be made shorter, which is 
definitely advantageous. Hence we search for pulses with 
lower amplitude among the second-order pulses. This is 
done by using an additional coefficient b„i in the ansatz 



( 33 1 and minimizing the amplitude Vq of the resulting 



solutions by varying this additional coefficient. 



A. First-order pulses 

For first-order pulses and the pure dephasing model, 
the set of conditions (21b I comprises only three equa- 
tions given in Appendix |A Adding conditions ( 34 1 and 



dephasing bath leading to six additional integral condi- 
tions besides the first-order terms (see Appendix|B]) . Note 
that more equations are to be fulfilled than for amplitude 
modulation 30J because the frequency-modulated pulses 
involve all three spin directions. The solutions for tt and 
tt/2 pulses are given in Table [Tl| and they are displayed in 
Figs. [3] and |4] Numerically, the double integrals in Eqs. 
( |B2[ ) are particularly demanding. Full quantum mechan- 
ical studies of higher-order pulses will be hampered by 
even higher dimensional integrals occurring in the Mag- 
nus expansion |38j . An alternative route, which may be 
numerically more efficient, consists of the direct solution 
of the Schrodinger equation |24j. The mathematical ex- 
istence of higher-order pulses is known "WT . 

Since we are interested in pulses with low amplitudes, 
we aim at minimizing the amplitude. To this end, we add 
another Fourier coefficient to the ansatz ( 33 1 and vary 



this additional parameter. In this way, we obtained the 
pulses FM-2-MIN-PI and FM-2-MIN-PI2 given in Table 
III and plotted in Figs. [5] and [6] Empirically it turned 
out to be more efficient to consider 614 instead of 611 as 
an additional coefficient. It is expected that even lower 
amplitudes can be achieved by using further coefficients. 
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Figure 6. (Color online) Minimized second-order n/2 pulse 
FM-2-MIN-PI2. We also plot sin$(i) oc Vy{t) and cos$(f) oc 
Vx{t) to illustrate the pulse shape in spin space. The left 
scale refers to ^{t) and the right scale to sin "I>(f) and cos ^{t), 
respectively. The Fourier coefficients for this pulse are given 
in Table Uni 



But our calculations with different coefficients, not sliown 
here, indicate tliat tliis route would improve tlie ampli- 
tude only by 1 — 2% at the expense of a more complex 
pulse shape. 

By using only one free coefficient (614) we found tt/2 
pulses with amplitudes lower than 9.0/rp to be compared 
with the amplitude-modulated pulses |30j with ampli- 
tude 11.5/Tp. For TT pulses we need Vq = 10.7/Tp in 
comparison to ll.O/rp for amplitude modulation. The 
amplitudes of the amplitude-modulated pulses refer to 
piecewise constant pulses; for continuous pulses they are 
even higher. Of course, the reduction of the amplitudes 
for frequency-modulated pulses is not spectacular. But 
it is remarkable that pulses with relatively low ampli- 
tudes can be found despite the larger number of con- 
ditions to be fulfilled: the frequency-modulated pulses 
fulfills 9 + 2 = 11 equations including the conditions for 
the angles, and the amplitude-modulated pulse fulfills 
5 + 1 = 6 conditions including the condition for the ro- 
tation angle. 



VI. CONCLUSIONS 

In this paper we extended the existing perturbative 
approach to decouple a spin from a quantum mechanical 
bath by means of short control pulses in two ways. 

First, we allowed for a time-dependent bath, which 
means both the bath Hamilton operator and the coupling 
operators may have an explicit, analytical time depen- 
dence. Yet, we found that this time dependence does not 
alter the requirements for the pulse shape which were de- 
rived previously for time- independent baths |30) . Hence, 
the pulses found previously are also applicable for time- 
dependent environments as they arise, for instance, in 



time-dependent reference frames or in the interaction pic- 
ture of otherwise time-independent Hamiltonians. This 
is our first key result. 

Second, we studied frequency-modulated pulses in first 
order and in second order in the pulse duration Tp for 
quantum mechanical baths. Previously, only amplitude 
modulation was considered explicitly for quantum me- 
chanical baths [TBI EH [30] . Frequency modulation was so 
far studied for static baths only [53] . We provide explicit 
solutions for continuous frequency-modulated pulses with 
amplitudes which have been minimized empirically. Such 
pulses are expected to be useful in experiments where 
no amplitude modulation can be realized or where the 
achievable accuracy for frequency modulation is supe- 
rior to the accuracy of amplitude modulation. For in- 
stance, they can be used to implement realistic optimized 
dynamic decoupling [321 IMl where the dynamic decou- 
pling sequence is adapted to pulses of finite length. The 
frequency-modulated pulses constitute our second key re- 
sult. 

We emphasize that modulated pulses correspond to 
quantum gates which are robust against decoherence in 
the framework of quantum-information processing, e.g., 
the 7r/2 pulse about ay preceded by a tt pulse about az 
realizes the important Hadamard gate up to a global fac- 
tor i. [iT] , 

Further work should concentrate on higher-order 
terms not studied here. Such terms comprise higher- 
dimensional integrals so that the numerical effort in- 
creases considerably. Another promising route is to ex- 
tend the model from pure dephasing to general decoher- 
ence. This would allow for systems with finite Ti as well, 
at the expense of more complex pulses. 

But at the present stage, it is also called for to verify 
the performance of the proposed pulses experimentally in 
order to assess how promising further extensions would 
be. 
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Appendix A: First-order conditions 



For the first-order conditions we inspect Eq. (21b I to 
find the corresponding scalar equations. In a purely de- 
phasing bath the sum over j collapses to j = z, resulting 
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in the scalar equations 



dephasing bath model and we eventually obtain 



Vii = / a,y{t) sin(-0(i)) - (1 - coa{'iJj))ax{t)az{t)dt 

(Ala) 

'712 = / a^it) sm{ip{t)) + (1 - cos{ip))ay{t)az{t)dt 



(Alb) 

m3= r cos(^(t)) + (1 - cos{ij))a,{tfdt. (Ale) 



The ai are the components of the global axis of rotation 
parametrized in Eq. ( 28 ) . 



and 



V2i= t [ay{t) sm{ipit)) - (1 - cos{ij))aa:{t)a^{t)] dt 
Jo 

(Bla) 

7722= / t[aj:{t)sm{^{t)) + {1 - cos{ij))ay{t)a;,{t)]dt 
Jo 

(Bib) 



??23 = / t 
10 



COS {ilj{t)) + (1 - cos{il}))a^{tY\ dt (Blc) 



7^24= / dti I dt2[ny^{ti)n^^{t2) - n22{ti)ny^{t2)] 
10 Jo 

(B2a) 

?725 = / dti / dt2[n^z{ti)nx^{t2) - n^z{ti)nzz{t2)] 
Jo Jo 

(B2b) 

l-rp r-ti 

ri2&^ I dii / dt2[nxz{ti)nyz{t2) - nyz{ti)n3:z{t2)] ■ 

(B2c) 



Jo 



The matrix elements nij{t) occurring here are those of 
the rotation matrix Dai—ip) given explicitly in Eq. (CI I. 
The components Ui are parametrized in Eq. (28|. 



Appendix B: Second-order conditions 



For the second-order conditions, we additionally have 



to consider Eqs. (22) to find the corresponding scalar 



equations. Again certain sums collapse due to the purely 



Appendix C: Rotation matrix 

To derive the matrix £)„(—'!/;) we refer the reader to 
Ref. (3U! It is calculated by comparison of the coefficients 



in Eq. (17a). We obtain the matrix (CI) below, where 
the time dependencies of ip(t) and d{t) are omitted for 
clarity: 



cos ^(J + {1 — cos i^)a'^ 



1 1/' + (1 — cos 'ip)axay —CLy sin + (1 — cos il>)o-x0.z 



Da{—''P)= —Qz siiiip + {1 ~ cos ip)axay cos V' + (1 — cos-0)a^ Ux simjj + (1 — cos ip)0'yO'z 
Qy s'mip + {1 — cos il;)axaz — Qx s'mip + {1 — cos ip)ayaz cos "0 + (1 — cos 7/^)0^ 



(CI) 
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